C inelastic scattering studied with microscopic coupled-channels method 
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In order to test the 16 C internal wave function, we perform microscopic coupled-channels (MCC) 
calculations of the 16 C(0f -» 2+) inelastic scattering by 208 Pb target at E/A=52.7 MeV using the 
antisymmetrized molecular dynamics (AMD) wave functions of 16 C, and compare the calculated 
differential cross sections with the measured ones. The MCC calculations with AMD wave functions 
reproduce the experimental data fairly well, although they slightly underestimate the magnitude of 
the cross sections. The absolute magnitude of calculated differential cross sections is found to be 
sensitive to the neutron excitation strength. We prove that the MCC method is a useful tool to 
connect the inelastic scattering data with the internal wave functions. 

PACS numbers: 24.10.Eq, 25.60.-t 



I. INTRODUCTION 

Recently, the opposite deformations between proton 
and neutron densities in C isotopes were theoretically 
suggested 0,0] by the method of antisymmetrized molec- 
ular dynamics (AMD) : The proton density has an oblate 
deformation, while the neutron density has a prolate de- 
formation, and the symmetry axis of proton is perpen- 
dicular to that of neutron in 10 C and 16 C. Based on 
this picture, the author also gave a qualitative expla- 
nation Q for unusually small electric transition strength 
B{E2-2+ -> Of) = 0.63 ± 0.19 e 2 fm 4 in 16 C, which is 
derived from life time measurement [4J. According to 
Ref.Jjjj], the 2^ state is a rotational excited state, and the 
rotational axis is perpendicular to the neutron symmetry 
axis. In this excitation mechanism, the proton transition 
strength is reduced due to the difference of deformation 
between the proton and neutron distributions mentioned 
above, and therefore, the 0^ — > 2^ transition is domi- 
nated by the neutron excitation. 

In order to search for the possible difference of proton 
and neutron contributions to excitation of the 2f state 
in 16 C, inelastic scattering experiment of 16 C on 208 Pb 
target was performed [5| applying the Coulomb-nuclear 
interference method. The analysis was carried out by us- 
ing the deformed potential model, and the proton and 
neutron transition matrix elements, M p and M n , were 
extracted. In Ref.p|, it is mentioned that the experimen- 
tal transition probability is inconsistent with theoretical 
ones (the AMD, extended AMD and shell- model calcula- 
tions are cited). However, the phenomenological analysis 
done in Rcf. I contains some assumptions, and hence, it 
seems inappropriate to compare the M p and M n values 
evaluated in Ref.[j| with those calculated theoretically. 

To test the 16 C internal wave function, we should di- 
rectly link the cross section with the wave function by 
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calculating the differential cross sections of the inelastic 
scattering of 16 C on 208 Pb target with the microscopic 
coupled-channels (MCC) method, and compare the cal- 
culated result with the experimental reaction data mea- 
sured in Ref.|5|. 

The MCC method has been applied for studying reac- 
tions of stable nuclei, such as the 6,7 Li and 9 Be elastic 
and inelastic scattering ||J Q, the resonance reactions 
of the 12 C+ 12 C system leading to inelastic j^, and 
8 Be+ 16 a-transferred channels |l0j, the rainbow scat- 
tering of 16 0+ 16 system ^jji e t c -; adopting the mi- 
croscopic cluster model wave functions of 6,7 Li 0, 9 Be 
[13, 12 C [3, and 16 14]. Because the microscopic 
cluster model wave functions well reproduce the mea- 
sured charge form factors not only of the elastic electron- 
scattering but also the inelastic one, the wave functions 
are reliable and adopted for studying nuclear reaction 
mechanisms. The MCC calculations successfully repro- 
duce the experimental reaction data. The reliability of 
the method has already been established. Hence, we 
think it is possible to examine inversely the validity of 
calculated internal wave functions by comparing the re- 
sult of MCC with experimental reaction data. 

In this paper, we adopt the AMD internal wave func- 
tion of 16 C. The reason why we think it is worthy to test 
the validity of the AMD wave function in the present 
study is as follows. (1) Because no inert cores and no 
clusters are assumed, the AMD wave function is flexible. 
Therefore, AMD is suited for structure study of general 
unstable nuclei. The applicability has been proved in 
many works |2(. (2) Deformations of proton and neu- 
tron densities are obtained dynamically. In other words, 
electromagnetic transition probability can be calculated 
without introducing effective charge. (3) In AMD, it is 
easy to carry out the spin-parity projection, which is nec- 
essary for microscopic calculation of transition density. 

It should be noted that any nuclear structure mod- 
els are applicable to MCC, if the diagonal density and 
transition ones to excited states can be calculated from 
a model wave function that gives no spurious center-of- 
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mass component. In the next section, we briefly describe 
the MCC method. More detailed description is made in 
Refs-jHH. 

II. FORMALISM 

A. Coupled-channels formalism 

The coupled-channels equation describing the collision 
of two nuclei for a total angular momentum of the system 
J is written as 

' T R + U<£ iaL (R)-E a ] xi J kR) 

= E u «l,«'l'(R) x%{R), (i) 

(a'L')^(aL) 

where Tr denotes the kinetic-energy operator, and a and 
L denote a channel and the orbital angular momentum 



associated to the relative coordinate R. In the present 
study, we take into account the elastic and 16 C excita- 
tion channels, while only the ground state (0 + ) is con- 
sidered for the target 208 Pb nucleus. Thus, the chan- 
nel a is designated by the spin I and the excitation 
energy e a of 16 C. E a represents the center-of-mass en- 
ergy for the projectile-target relative motion in channel 

a (E a = .Ec.m. — £<*)■ xUi'(-R) is the relative wave func- 
tion and is obtained by solving Eq.Q numerically. 



In Eq.Q, V a l a , L , (R) represents the diagonal (a,L) = 
(a', L') or the coupling (a, L) ^ (a', L 1 ) potential, which 
is composed of the nuclear part V^ J ^, L ,(R) and the 

C( J) 

Coulomb part V a ^ J, L , (R). The nuclear part is given 
by the double-folding model and defined as 



V£U«) = ^E^'^ W(ILI'L';JX) (L0L'0|A0) 

A 

x L fdk dv.dv, [ Voo (x) (p?W(n) + fi^inj) (p^ 0) (r 2 ) + p^\r 2 )) 
(x = n+R-r 2 ) 



r 



(2) 



where W(ILI'L'; JA) represents the ordinary Racah co- 
efficient and I is y/2I + 1. p^j,(ri) and p™\, (r*i) 

G°oo^( r 2) and Poo ( r 2)) are the radial components of 
proton and neutron transition densities of 16 C ( 208 Pb), 
respectively, which will be mentioned in the next subsec- 
tion in detail. Doo( x ) represents the spin- and isospin- 
scalar (S = T = 0) component of an effective nucleon- 
nucleon interaction, while i>oi(x) represents the spin- 
scalar, isospin-vector (S — 0, T = 1) component. For 
this effective interaction, we adopt the DDM3Y (density- 
dependent Michigan three-range Yukawa) [H which 
is defined by 



■yoo(oi) (E, p; r) = 500(01) (E, r)f(E, p), 



(3) 



where r is the intcrnucleon separation, and f(E 1 p) is a 
density dependent factor 



f(E,p) = C(E) [l + a(E) e-PW<>] 



(4) 



Here, E denotes an incident energy per nucleon in the 
laboratory system. The coefficients C(E), a(E), and 
(3(E) in the density-dependent factor f(E,p) were de- 
termined at each energy by fitting a volume integral of 



the Voo (E, p; r) to the real part of the optical potential 
felt by a nucleon in the nuclear matter Il7 j . q pocn ii (E, r) 
in Eq.@ is the original M3Y interaction flSt tl9j|: 



-4r e~ 2 - 5r 

300 (E,r) = 7999- 2134—— 

4r 2.5r 



■Joo(E)6(t) MeV, 
(5) 



with 



and 



J 00 (E) = - 276 (1 - 0.005.E) MeV fm 3 , (6) 



— 4:T 2 5v 

goi(E,r) = -4886^— +1176^-— - 
4r 2.5r 



with 



■J i(E)6(r) MeV, 
(7) 



J m {E) = 228.4 (1 - 0.005S) MeV fm 3 . (8) 



The units for E and r are MeV/nucleon and fm, respec- 
tively. 

The Coulomb part V^ J J, L ,(R) is also given by the 
double-folding model. The double-folded Coulomb po- 
tential is written in the same form as Eq.(J2J), by replacing 
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TABLE I: B{E2), M { v i] and values of 16 C calculated 

by AMD, in which the strength of the spin-orbit force is set 
to (i) u Js =900 MeV, (ii) w is =1500 MeV, and (iii) u Js =2000 
MeV. The experimental data of B{E2) is taken from Ref.Q . 







(0 


(ii) 


(iii) exp. 


B(E2;2+ 


- 0+)(e 2 fm 4 ) 


1.9 


1.4 


0.9 0.63 ±0.19 


A4 2) (2+ 


- 0+) (fm 2 ) 


3.1 


2.6 


2.2 


Mi 2) (2+ 


- 0+) (fm 2 ) 


13.0 


12.2 


8.9 



the neutron densities and the nucleon-nucleon interaction 
as 

e 2 

«oo(x) -> — , 
uoi(x) -> 0. 

Since DDM3Y has no imaginary part, we add the imag- 
inary potential W^ J J, L ,(R) to the nuclear part, which is 

assumed as W^, L ,(R) = Nj ■ V^ L [ J J, L ,(R), where Nj 
is the only a phenomcnological parameter of the present 
MCC formalism. The simple assumption for the imag- 
inary part should be valid in the present case, since we 
only discuss the cross sections at very forward scattering 
angles, which is not sensitive to the detail shape of the 
potential in whole radial range. Hence, the interaction 
potential has the form as 

uS;%{R) = (1 + iNi)vj*$ v {R) + ( 9 ) 



B. Transition density 

The diagonal or transition density of proton at a posi- 
tion r with respect to the center-of-mass of the nucleus 
can be expanded into multipole components : 

z 

i, v W = (^(6 IE *( r - r 0li&„'(0) 

i=l 

= Y,(l'v'^\Iv) fi { v\r)Y^(r), (10) 

X.j-i 

where V'? l /(£) represents the proton wave function in the 
nucleus, fip (r) represents the radial component of the 
transition density, which is used in Eq.Q. The radial 

component of the neutron transition density p^j, (r) is 
obtained in the same manner as the proton case in terms 
of the neutron wave function The proton or neu- 

tron matrix element of rank A is defined as 

M W(/' -+I)=lJ p$\r) r x+2 dr, (11) 



where r represents p or n. The proton matrix element is 
related with the electric transition strength B(EX) as 

B (w^)= |M " ( ;r J)|2 - 2 < i2 > 

Here, we use the AMD wave function for tpj u (£,) and 
if>f v (C) to calculate the transition densities defined in 
Ea.lfTUjl. We consider the ground state (Of) and first 
excited 2f state. In Ref.Ql, two versions of 16 C internal 
wave function are obtained in the variation before projec- 
tion (VBP) formalism changing the strength of the spin- 
orbit force, (i) u/ s =900 MeV and (ii) w is = 1500 MeV. The 

(2) (2) 

B(E2) values as well as the M p ' and Mn ' ones obtained 
with AMD wave functions (i) and (ii) are summarized in 
Table U with the experimental data of B(E2). With the 
spin-orbit forces (i) and (ii), the AMD calculation repro- 
duces well the systematic behavior of the B(E2) value 
and root-mean-square radius of the C isotopes as shown 
in Ref.0- In particular, the systematic feature that the 
B(E2) value of 16 C is abnormally small compared with 
other C isotopes ( 10 C, 12 C and 14 C) is well reproduced 
by AMD with the spin-orbit (i) and (ii), although the 
B(E2) value of 16 C is slightly overestimated. In addition 
to above two, we also use the AMD wave function for 
which the strength of the spin-orbit force is set to (iii) 
m; s =2000 MeV so as to reduce the B(E2) value and to be 
close to the experimental value. The B(E2) value as well 

(2) (2) 

as the Mp ' and M}?> ones of the case (iii) is also shown 
in Table U One might think that the case (iii) gives the 
best wave function. However, increasing the strength of 
the spin-orbit force to reduce the B(E2) value of only 
16 C as done in the case (iii) may lead to an unrealistic 
situation, because the systematic behavior of the other 
C isotopes is not reproduced with such strong spin-orbit 
force. Hereafter, we refer the three versions of AMD wave 
function as AMD(i), AMD(ii) and AMD (iii). 

Figure shows the radial components of the diago- 
nal and transition densities obtained with AMD(i). In 
the upper panel (A), the A=0 components of the diag- 
onal density are shown. The solid curves represent the 
proton and neutron densities of the ground state (0^), 
and the dashed curves represent those of the 2f state. 
It is found that the shape of the 2f diagonal densities 
are almost the same as those of Of , except for the region 
around the origin. The proton density of the 2f state 
is almost identical to that of the 0^ state and the dif- 
ference cannot be seen in the figure. In the lower panel 
(B), the A=2 components of the diagonal and transition 
density are shown. The dotted and solid curves represent 
the proton and neutron transition densities, respectively, 
for the 2~\ — * Of transition. The transition is found to 
be dominated by the neutron component especially in 
the surface region. The double-dot-dashed and dashed 
curves are the A=2 components of the 2f diagonal den- 
sity for the proton and neutron. These two curves are 
displayed with opposite sign in the figure. The shapes 
of the double-dot-dashed and dashed curves are similar 
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FIG. 1: The radial components of the diagonal and transition 
densities obtained with AMD(i) for 16 C. (A) A=0 component. 
The solid curves represent the proton and neutron diagonal 
density of the ground state (Of), while the dashed curves 
represent the those of the state. The proton density of 
the 2^ state is almost identical to that of the OJ" state and 
the difference cannot be seen in the figure. (B) A=2 compo- 
nent. The dotted and solid curves represent the proton and 
neutron transition densities from 2^~ to 0^~. The double-dot- 
dashed and dashed curves represent the X—2 components of 
the proton and neutron diagonal densities for the 2J" state. 
These two curves are displayed with opposite sign. 



to the dotted and solid curves, respectively. It should be 
noted that the proton part of X—2 component of the di- 
agonal density is proportional to the electric quadrupole 
moment of the 2+ state. We neglect the A=4 component 
of the 2f diagonal density, because this component is 
very small and is expected to have a small contribution 
to the inelastic scattering. 

In Fig|2lA), the proton transition density of AMD(i) 
(solid curve) is compared with those of AMD(ii) and 
AMD(iii), which are represented by the dotted and 
dashed curves, respectively. The transition densities of 
AMD(ii) and AMD(iii) are smaller than that of AMD(i), 
and the transition density of AMD(iii) is slightly shifted 
to small r side. The difference of the proton transition 
density causes the difference of the electric transition 
strength, as shown in Table [I] In Fig^B), the neutron 
transition densities of AMD(i), (ii) and (iii) are shown 
by the solid, dotted and dashed curves, respectively. Al- 



FIG. 2: Comparison of the proton (A) and neutron (B) tran- 
sition densities for 2^ — > 0^. The solid, dotted and dashed 
curves represent the transition densities of AMD(i), AMD(ii) 
and AMD (iii), respectively. The solid curves in (A) and (B) 
are the same as the dotted and solid curves in Fig^B), re- 
spectively. 



though the overall shapes of transition densities of (i) and 
(ii) are almost the same, the magnitudes are found to be 
slightly different. The magnitude of transition density of 
AMD (iii) is suppressed when it is compared with those 
of (i) and (ii), especially in inner region. 

The proton density distribution of the 208 Pb nucleus is 
obtained by unfolding the charge density 20] , which was 
obtained by a model-independent analysis of electron- 
scattering experiment, with the realistic proton charge 
form factor [21|. The neutron density distribution is ob- 
tained by assuming that the shape of the neutron den- 
sity is the same as that of proton one, namely p n (r) = 
(N/Z)p p (r). This assumption is known to be valid for 
stable nuclei. 

In the next section, we show the results of the MCC 
calculation using the AMD transition densities described 
above. 



III. RESULTS 

In Fig|31 we show the double-folding model potential 
using the densities obtained by AMD(i) described in sub- 
section ^TOJ The solid curve represents the nuclear po- 
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FIG. 3: The diagonal (A) and coupling (B) potentials of the 
16 C+ 208 Pb system calculated with AMD(i) are represented 
by the solid curves. The dotted and dashed curves are the 
contributions of the proton and neutron components of 16 C, 
respectively, where both the proton and neutron components 
of 208 Pb(0+) are included. 



tential Eq.JSJ. The dotted and dashed curves are the 
contributions of proton and neutron components of 16 C, 
respectively, where both the proton and neutron com- 
ponents of 208 Pb(0^~) are included. In FigEjA), the 
diagonal potential of 16 C(0^) + 208 Pb elastic channel 
is shown. Since the diagonal density of the 2\ state 
in 16 C resembles closely that of Of state as shown in 
FigfflA), the diagonal potential of 16 C(2+) + 208 Pb is 
almost the same as that of 16 C(0| ) + 208 Pb, and there- 
fore, not shown here. In Figl^B), the coupling potential 
of 16 C(0^ -> 2+)+ 208 Pb is shown. It is found that the 
neutron component has a dominant contribution to the 
total potential, especially in vicinity of the strong absorp- 
tion radius rsA ~ 11 fm. 

In order to see the effect of the isovector component, 
we also decompose the diagonal and coupling potentials 
into the isoscalar and isovector components. The result 
is shown in FigQ] The solid curves are the same as the 
solid ones in FigEl and the dotted and dashed curves are 
the isoscalar (IS) and isovector (IV) components, respec- 
tively. Compared with the isoscalar component, the mag- 
nitude of the isovector component is about 5%, and the 
sign is opposite. The 5% reduction of the coupling poten- 
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FIG. 4: The isoscalar (dotted) and isovector (dashed) com- 
ponents of the diagonal (A) and coupling (B) potentials. The 
solid curves are the same as those in FigE] 



tial leads to about 10% reduction of the inelastic scatter- 
ing cross sections and is not negligible in the present case. 
Therefore, we include the isovector component through- 
out the present calculations. 

We perform a coupled-channels calculation using the 
diagonal and coupling potentials shown in FigEI Because 
of the high incident energy, the coupled-channels equa- 
tions are solved numerically with the relativistic kine- 
matics, which has a non-negligible effect at forward-angle 
cross sections. The parameter for the imaginary poten- 
tial Ni is set to 1.2. Following the procedure of Ref.pJ, 
the calculated cross sections are smoothed by Gaussian 
functions according to the experimental angular uncer- 
tainty of 0.28°. The result is shown in FigGO The differ- 
ential cross sections are shown as a function of scattering 
angle 9 in the laboratory system. The crosses are the 
experimental data and the solid curve represents the 
result of the coupled-channels calculation. It is found 
that the MCC calculation with AMD(i) reproduces the 
experimental data fairly well, although it slightly under- 
estimates the magnitude of the cross sections at large 
angles. While the oscillatory shape of angular distribu- 
tion is formed by interference between the nuclear and 
Coulomb components, which are represented by the dot- 
ted and dashed curves, respectively, in FigG3 the average 
strength of the calculated cross section is determined by 
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FIG. 5: Angular distribution of 16 C(0^ -> 2+) inelastic scat- 
tering on 208 Pb target at E/A=52.7 MeV. The crosses are the 
experimental data and the solid curve is the result of our cal- 
culation. The dotted and dashed curves are the nuclear and 
Coulomb components, respectively. 
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FIG. 7: (A) The coupling potentials calculated with AMD (i) , 
(ii) and (iii) are shown by the solid, dotted and dashed curves, 
respectively. (B) The results of the MCC calculations with 
AMD(i), (ii) and (iii) are shown by the solid, dotted and 
dashed curves, respectively. The crosses are the experimental 
data. 
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FIG. 6: Ni dependence of the calculated differential cross 
sections. The dotted, dashed, solid, dot-dashed and double- 
dot-dashed curves are the results of the MCC calculations 
with Ni — 1.0, 1.1, 1.2, 1.3 and 1.4, respectively. 



the nuclear excitation. Particularly, the neutron com- 
ponent dominates the nuclear excitation as understood 
from FigEJtB). Therefore, present result indicates that 
AMD(i) slightly underestimates the neutron excitation 
strength by about 10 %, while it overestimates the pro- 
ton excitation one as shown in Table [fl 

The parameter iVj cannot be determined theoretically 
in the present MCC framework. In order to see the Nj 
dependence of the calculated result, we perform the same 
calculation with different Nj values. The results are 
shown in FigEJ The dotted, dashed, solid, dot-dashed 
and double-dot-dashed curves are the results of the MCC 
calculations with iVj=1.0, 1.1, 1.2, 1.3 and 1.4, respec- 



tively. It is seen that A^j-dependence is very weak, al- 
though the cross sections at very forward angles slightly 
changes with Nj. The angular distribution around 6>i a b 
= 3-4 degrees is seen to be independent of the Ni 
value. Because the calculation with Nj— 1.2 reproduces 
the shape of the data a little better than others, we 
choose JVj=1.2 in the present calculation. 

Next, we perform the MCC calculations using AMD(ii) 
and (iii). In Fig|J{A), the coupling potential with 
AMD(i) represented by the solid curve is compared with 
those with (ii) and (hi), which are represented by the 
dotted and dashed curves, respectively. It is noticed that 
the strength of the coupling potential is almost propor- 
tional to the M^ 2) + Mi 2) value; 16.1 fm 2 for (i), 14.8 
fm 2 for (ii) and 11.1 fm 2 for (iii). The results of the 
MCC calculations with AMD(ii) and (iii) are shown in 
FigJJfB) by the dotted and dashed curves, compared with 
the result of AMD(i) shown by the solid curve which is 
the same as that in Fig|Sl and the experimental data. 
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As expected from the strength of the coupling potential, 
the differential cross sections with AMD(ii) are slightly 
smaller than those calculated with AMD(i), and those 
with AMD(iii) severely underestimates the magnitude of 
the measured cross sections. The magnitude of differen- 
tial cross section of inelastic scattering directly reflects 
the electric and hadronic transition strength of the 16 C 
nucleus. For AMD(iii), the proton transition seems good 
because it gives the B(E2) value being close to the mea- 
sured one. However, the nuclear excitation strength of 
(iii) is too small as shown in Fig[7fB). Since the nuclear 
excitation is dominated by the neutron component, this 
result indicates that AMD (iii) fails to describe the neu- 
tron excitation correctly. This fact cannot be known by 
the experimental data of electromagnetic probe. There- 
fore, we think it is very important that the internal wave 
function of a nucleus obtained theoretically is tested by 
the experimental data of hadronic probe to investigate 
the behavior of neutron component. 

For further investigation, it is very interesting to ana- 
lyze the experimental data of the 16 C inelastic scattering 
on proton measured at RIKEN , which is more sensi- 
tive to the neutron excitation. 

Finally, we see how the angular distribution changes if 
the symmetry axis of proton is aligned to that of neutron. 
In the original AMD wave function, the symmetry axis 
of proton is perpendicular to that of neutron. We arti- 
ficially make the aligned proton density from the AMD 
wave function (ii) by rotating the proton wave function 
to set its symmetry axis to be parallel to the neutron 
one. The aligned proton density is shown in FiglSJA) by 
the dotted curve. The solid curve is the neutron tran- 
sition density which is the same as the dotted curve in 
Fig|2lB). The proton transition density has the opposite 
sign to the neutron one in the surface region. The result 
of the MCC calculation using this aligned proton tran- 
sition density is shown in Fig|SjB). It is found that the 
calculated angular distribution is out-of-phase compared 
with the experimental data. This result indicates that 
the proton transition density should have the same sign 
as the neutron transition density in the surface region. 



IV. SUMMARY AND CONCLUSION 

In order to test the 16 C internal wave function, we 
studied the 16 C(0^ 2^)inelastic scattering on 208 Pb 
target at E/A=52.7 MeV |f| by the microscopic coupled- 
channels (MCC) method using the internal wave func- 
tion of the 16 C nucleus obtained by the antisymmctrizcd 
molecular dynamics (AMD) 3]. In Ref.0, two versions 
of wave function are obtained with the strength of spin- 
orbit force (i) u ls = 900 MeV and (ii) u&=1500 MeV. 
It was shown in Ref.0] that these AMD calculations re- 
produced the systematic behavior of B(E2) value and 
root-mean-square (RMS) radius of C isotopes. The MCC 
calculations using these wave functions of 16 C reproduce 
well the measured differential cross sections, although 
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FIG. 8: (A) The proton transition density when the symmetry 
axis is artificially rotated to be parallel to the neutron one is 
shown by the dotted curve. The neutron transition density 
shown by the solid curve is the same as the dotted curve in 
FigUlB). (B) The result of the MCC calculation when the 
aligned proton density shown in (A) is used. 



they slightly underestimate the magnitude of the cross 
sections at large angles. Especially, the shape is rather 
well reproduced around 0i a t>= 3-4 degrees, where the an- 
gular distribution is independent of the strength parame- 
ter Ni of the imaginary potential. While the shape of dif- 
ferential cross section due to the interference between the 
nuclear and Coulomb excitation components is sensitive 
to the strength of the proton excitation, the magnitude 
of the cross section is sensitive to the strength of neutron 
excitation, because the nuclear excitation is dominated 
by the neutron one in the present case. Therefore, we 
can conclude that the AMD wave function of (i) predicts 
the neutron excitation strength of 16 C reasonably well, 
although the strength may be slightly underestimated by 
about 10 %. 

We also performed a coupled-channels calculation us- 
ing AMD wave function, for which the strength of the 
spin-orbit force is set to (iii) 1(^=2000 MeV. Although 
this wave function gives B(E2) value of 16 C close to the 
measured one, the systematic behavior of B(E2) and 
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RMS radius of the other C isotopes are failed to be re- 
produced due to the unrealistic strength of the spin-orbit 
force. The MCC calculation using AMD(iii) severely 
underestimates the differential cross sections, which in- 
dicates that the neutron excitation is not properly de- 
scribed when the spin-orbit force (iii) is used. It can 
be said that testing the validity of calculated wave func- 
tion only with the electromagnetic experimental data, 
such as B(E2) value, may be insufficient. Especially for 
neutron-rich nuclei, it is expected that the proton density 
is largely different from the neutron one. Therefore, we 
think it is very important that the internal wave func- 
tions of a nucleus obtained by any nuclear structure the- 
ory should be tested not only by the experimental data of 
electromagnetic probe but also those of hadronic probe, 
particularly to investigate the behavior of the neutron 
component, as done in the present paper. 

We showed that the MCC calculation is a useful tool to 
link the inelastic scattering data with the internal wave 
functions obtained theoretically. Note that since the di- 
agonal density is also reflected by the behavior of the 



calculated differential cross section through the diago- 
nal potential as the distorting effect, the overall feature 
of diagonal and transition densities of nucleus can be 
tested in the consistent procedure. Nuclear reaction data 
themselves are available for the nuclear structure study 
similarly to the RMS radius, electromagnetic transition 
strength, charge form factor, etc., by applying the MCC 
method. 
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